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We study coherent photoassociation, phenomena analogous to coherent optical transients in few- 
level systems, which may take place in photoassociation of an atomic Bose-Einstein condensate 
but not in a nondegenerate gas. We develop a second-quantized Hamiltonian to describe photoas- 
sociation, and apply the Hamiltonian both in the momentum representation and in the position 
representation (field theory). Solution of the two-mode problem including only one mode each for 
the atomic and molecular condensates displays analogs of Rabi oscillations and rapid adiabatic pas- 
sage. A classical version of the field theory for atoms and molecules is used to demonstrate that, in 
the presence of photoassociating light, a joint-atom molecule is unstable against growth of density 
fluctuations. Experimental complications, including spontaneous emission and unwanted "rogue" 
photodissociation from a photoassociated molecule are analyzed. A two-color Raman scheme is 
studied as a method to set up an effective two-mode scheme with reduced spontaneous emission 
losses. We discuss photoassociation rates and photoassociation Rabi frequencies for high-lying vi- 
brational states in alkalis both on the basis of molecular-structure calculations, and by comparing 
with an experiment [Wynar et al, Science 287, 1016 (2000)]. 
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I. INTRODUCTION 

Laser cooling and its spin-offs now routinely pro- 
duce gaseous samples in which thermal energies, when 
expressed as frequencies, are smaller than the typical 
linewidth of an optical dipole transition. Photoassoci- 
ating transitions, in which two thermal atoms combine 
in the presence of light to make a molecule, may there- 
fore exhibit linewidths every bit as narrow as the tran- 
sitions one encounters in nonlinear laser spectroscopy. 
As a result, photoassociation spectroscopy 0] has be- 
come the source of the most accurate molecular struc- 
ture data available. Bose-Einstein condensation j| is 
another recent triumph in the quest toward low temper- 
atures in atomic physics. The connection between pho- 
toassociation and Bose-Einstein condensation has long 
been close, though somewhat incidental; photoassocia- 
tion spectroscopy has provided key numerical data for 
condensation experiments ||. 

There have been early discussions of photoassociation 
of a condensate itself jl| ||. However, to us, the true 
scope of the connection between photoassociation and 
condensation was only revealed by our explicit observa- 
tion H that in a thermal gas it is the same phase space 
density that governs both the onset of Bose-Einstein con- 
densation and the efficiency of photoassociation. 

At the heart of our photoassociation work lies the qua- 
sicontinuum (QC) approach || |7|]. The idea is to enclose 
two colliding atoms in a box, which has the effect of dis- 
cretizing the dissociation continuum of the correspond- 
ing diatomic molecule. At the end of the calculations, 
the quantization volume is taken to infinity. Aside from 
resolving certain mathematical difficulties, this method 
turns out to have the unexpected benefit that analysis of 
photoassociation is reverted to analysis of few-level sys- 
tems, as in quantum optics or laser spectroscopy. For 
instance, studies of two-color photoassociation schemes 
may draw from decades of experience in quantum optics 



and laser spectroscopy (jfl. 

However, in its initial form our QC approach does not 
apply to a quantum degenerate sample. We sought to rec- 
tify this shortcoming by introducing a phenomenological 
second-quantized Hamiltonian for photoassociation [0. 
This idea was developed at the same time independently 
by Drummond et al. 0, and mathematically closely re- 
lated approaches to the Feshbach resonance are also un- 
der active study jlC| O, [l2|. Comparison with the QC 
approach gives the transition matrix elements to insert 
into our Hamiltonian. 

We first considered a two-mode model that only takes 
into account one CM. wave function for atoms and 
one for molecules, the modes containing the atomic and 
molecular condensates The main finding was coher- 
ent photoassociation analogous to coherent transients in 
few-level systems. For instance, the system may exhibit 
a form of Rabi flopping between atoms and molecules. 
Moreover, by properly sweeping the frequency of the 
photoassociating laser, in a process akin to rapid adia- 
batic passage the atomic condensate may be turned into 
a molecular condensate ||. In another development in 
this direction, we have argued that two-color free-bound- 
bound stimulated Raman adiabatic passage, STIRAP, 
is feasible starting from an atomic condensate |0|. We 
have also gone beyond the two- and three-mode approx- 
imations, allowing for an arbitrary position dependence 
of the atomic and molecular condensates, albeit in a clas- 
sical approximation similar to the one underlying the 
Gross-Pitaevskii equation iQ. R then turns out that 
an equilibrium with both atomic and molecular conden- 
sates present together with the photoassociating light is 
unstable. The sample tends to collapse spontaneously 
into clumps whose densities increase with time JT^ |. 

The primary purposes of the present paper are to 
document the numerous technical and physical details 
of our second-quantized approach to photoassociation 
that could not be accommodated by the letter format 
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of Refs. H , Jl3| , and |Q , and to extend our discussion 
in several directions that support those references. To 
offer a comprehensive account of the field theory version 
of our approach to coherent photoassociation, we have 
found it necessary to analyze the dipole matrix clement 
for photoassociation in detail. This endeavor in effect 
constitutes an alternative derivation (c.f. Ref. Q) of 
our entire QC methodology. Second, we add an analysis 
of two-color photoassociation of a quantum-degenerate 
sample in a three-mode approximation that is to some 
extent complementary to the one in Ref. |13[ . Much 
as expected, the two-color scheme provides a reprieve 
from spontaneous-emission losses from the primary pho- 
toassociated state. Third, we present a quantitative 
analysis of "rogue" photodissociation from a molecu- 
lar condensate to atomic modes outside the condensate. 
Our suggestion || that with increasing light intensity 
the unwanted photodissociation may overtake coherent 
condensate-condensate transitions is corroborated. We 
find a minimum usable time scale proportional to the 
the inverse of the recoil frequency of laser cooling. 



II. FIELD THEORY FOR ATOMS AND 
MOLECULES 

The task of the present section is to develop in detail 
the second-quantized approach governing photoassocia- 
tion of atoms into molecules in the prototype case of one 
laser color only. Simple heuristic arguments based on 
Refs. 10, [?j could, and did, achieve most of our aims in 
Ref. B where we dealt primarily with the momentum 
representation, but the field theory of Ref. |l4| calls for a 
few additional angles. To support them, we present here 
a partially new ab-initio discussion of our QC method. 

We take the photoassociating atoms to be in precisely 
one internal state, and similarly we assume that photoas- 
sociation leads to molecules with precisely one internal 
state. These assumptions could be relaxed, but then one 
has to follow the fate of the internal states as well. We do 
not go into this, but in essence assume that (i) the atoms 
are polarized and that (ii) the photoassociation resonance 
in itself selects a unique final state for the molecule. 



Probably the most prominent qualitative finding 
emerging from our work is the observation that it is Bose 
enhancement that ultimately facilitates coherent tran- 
sients such as Rabi flopping, adiabatic following, and 
ST1RAP in photoassociation of a condensate [^3[ fl5| . 
Throughout this paper, we continue to demonstrate how 
coherent optical transients come about in a condensate, 
and argue why they should be absent in a nondegenerate 
gas. 

There has recently been a remarkable experiment on 
two-color photoassociation of a condensate |r| . Accord- 
ingly, we include a detailed discussions on the values of 
experimental parameters in alkalis in general, and the pa- 
rameters of Ref. in particular. The analysis of the 
actual experiment demonstrates that there is still some 
way to go before genuinely coherent photoassociation is 
reached. 



In Sec. II we give a walk-through of our second- 
quantized Hamiltonian, including a detailed discussion of 
the dipole moment matrix element and both the momen- 
tum and position representation of the Hamiltonian. The 
special case with only one spatial mode for both atoms 
and molecules is the subject of Sec. [II. The classical ver- 



sion of the field theory, but including all spatial modes, 
is the subject of Sec. |y|. The numerous complications 
to our one color scheme that one is liable to encounter in 
real experiments, as well as the experimental parameter 
values, are the subject of Sec. |v|. The brief remarks in 
Sec. conclude the paper. There are also two appen- 
dices70 on the details of the relation of the dipole matrix 
element between second-quantized and quasicontinuum 
approaches, and [b| on the role of atom-atom collisions in 
our development. 



A. Two atoms 

We begin with a pair of atoms, assumedly in the dis- 
sociation continuum of a given potential energy curve of 
a diatomic molecule. As the atoms interact, their rela- 
tive momentum need not be a constant of the motion. 
Nonetheless, given a finite range for atom-atom interac- 
tions, in free space the wave functions of the relative mo- 
tion 0k (r) could still be characterized by the asymptotic 
(r — > oo ) wave vector k. 

On the other hand, in the spirit of the QC method j^, 
0, we assume that the relative motion of the atoms is 
confined to a finite volume V. There are two basic ques- 
tions in our two-atom analysis that must be considered. 
First, our phenomcnological many-particle Hamiltonian 
(Ref. M and Eq. ( |l2| ) below) is written down in terms of 
plane waves, yet it is more common to analyze photoas- 
sociation in terms of angular-momentum partial waves. 
How do we make the connection? Second, in our qua- 
sicontinuum method we resort to a finite quantization 
volume V , which tends to infinity only at the end of the 
calculations. How should we handle the finite quantiza- 
tion volume? 

We quantize the relative motion of the two atoms in a 
spherical box of radius R and volume V — ^R 3 using re- 
flecting boundary conditions. While angular momentum 
is still a constant of the motion, the usual eigenstates Q| 
of the relative motion cannot be characterized by a mo- 
mentum vector, even asymptotically. Nonetheless, it is 
evidently possible to construct orthonormal superposi- 
tions 0k (r) of the eigenstates of the spherical box that in 
the limit of a large box turn into the states 0k (r), i-e., 
states that behave like plane waves at large distances. 

Let us first take the inner product of a true plane wave, 
normalized to the volume V, with a spherically symmet- 
ric (real) test function f(r) with a finite range <C R. In 
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the limit k^Owe have 



ii = lim — = / d A r e l 
k-o y/V Jo 



7(r)~-=/ 47rrUrf(r) 
vV Jo 



(1) 

Next we take the same inner product for the £ = partial 
wave of the plane wave e lkr , basically the spherical Bessel 
function jo(kr) cx sin kr/r. We normalize this partial 
wave also in the spherical box, or equivalently, in the 
radial coordinate r with respect to the measure Airr 2 dr. 
In the limit of small k we have the integral for the inner 
product 
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k 



Obviously, the ratio of the two, 
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is the expansion coefficient of the s-wave in the plane 
wave, given that both the plane wave and the s-wave are 
normalized to the volume V. Actually, the k values of 
the spherical eigenmodes are quantized, and the smallest 
value is k = ir/R. We therefore have to be careful with 
the limit k — > when applying Eq. (^). 

When atom- atom interactions are taken into account, 
the radial eigenstates for a given angular momentum are 
not just spherical Bessel functions, but also reflect the 
scattering phase shifts. Henceforth we assume that only 
s-wave scattering needs to be considered, as is usually 
the case for bosons at sufficiently low temperatures. Ob- 
viously, even though the partial waves have phase shifts, 
by tweaking the phases of the partial waves it is_still pos- 
sible to make an asymptotic near-plane wave 0fc(r) out 
of the eigenstates of the relative motion of the atoms in 
the sphere. The adjustment of phases has no effect on 
the weight of the s-wave component in 4>k (r) , which may 
still be inferred from Eq. (j^). 

Let us denote the wave function of the particular 
molecular state we are aiming for by ip(r), and make the 
standard (albeit crude) approximation that the relevant 
electronic dipole matrix element d is a constant indepen- 
dent of the relative coordinate r of the atoms comprising 
the molecule. Then the QC dipole matrix element for 
photoassociation and photodissociation characterizing a 
process in which the relative momentum of the colliding 
atoms is Kk reads simply 



d(k) =d / d\ ^*(r)0 k (r). 



(4) 



We only consider s-wave collisions, and correspondingly 
set ip( r ) = "0( r ) as appropriate for a nonrotating J = 
molecule. Of course, only the / = component of 



the wave function <i>k(r) counts. Taking into account the 
weight from Eq. (|3|), we have 



47rr 2 drV*(r)^(r). (5) 



The notation <fP k (r) stands for the radial I = wave func- 
tion of the relative motion of the two atoms correspond- 
ing to the wave number k, and normalized to volume V 
as usual. 

A convenient qualitative model is provided by the lim- 
iting form 



4>%{r) = 



1 sin k(r — a) 



(6) 



where a is the s-wave scattering length. Combination of 
(||) and (H) gives 



d(k) 
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r dr sin k(r — a)ip* (r) . (7) 



The form (g) is only valid outside the range of the atom- 
atom interaction potential, and for small enough k so 
that the scattering phase shift may be written as ka. 
Equation (Q) is therefore quantitatively reliable only if 
the vibrational bound-state wave function of the molecule 
V>(r) happens to reside outside the range of the atom- 
atom interactions of the photoassociating atoms. In fact, 
this is the case in some of the current photoassociation 
experiments 17, fl8| . 

The matrix element (R) is explicitly proportional to 
V~ 1//2 . In fact, the V~~ 2 scaling of d(k) is a generic 
property of our QC approach, and does not depend on 
the specific assumptions of Eq. (Q). First, the almost- 
plane wave 4>k(r) is normalized to the volume V. No 
matter what kind of (square integrable) structure occurs 
around r ~ 0, the plane- wave form at large distances 
makes the normalization constant proportional to V~ x l 2 . 
Second, because of the finite range of the bounded molec- 
ular wave function ip( r )i its normalization coefficient does 
not depend on V. Third, thanks to the same finite range, 
the integral (Q) effectively extends only over a finite, fixed 
volume. The net result is that the matrix element indeed 
is proportional to the normalization constant of </>k(r), 
ocV- 1 / 2 . 

Another important property of the matrix element (Q), 
that it tends to a constant as k — > 0, is also a generic 
feature of the physics. Namely, for small enough k, the 
form of the s-wave (p^if) indeed is sinfc(r — a)/r, except 
within the range of atom-atom interactions. But within 
this range, the shape of (f>1(r) is independent of A; at small 
enough k. In order that the inner and the outer forms join 
smoothly, the amplitude of the inner wave function must 
be oc k for small k, just as is the amplitude of the outer 
wave function. So, for small enough k, the wave function 
</>°(r) is oc k for all of those r for which the bound-state 
wave function ip(r) is effectively nonzero. The matrix 
elements (|J), (|^) and (Q) therefore all tend to a constant 
as k — > 0. 
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In order to prepare for the comparison between pho- 
toassociation and photodissociation, let us assume that 
a laser field with the amplitude 



E(t) 



(8) 



is incident on a molecule. We denote the detuning of the 
light above the photodissociation threshold by S. Given 
the reduced mass of the colliding atoms fi, the corre- 
sponding resonant wave number kg and velocity i>o are 
such that /zfep /2/i = S and vq — hko/n- Using the stan- 
dard dipole and rotating-wave approximations as well as 
Fermi's golden rule, we have the photodissociation rate 
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(10) 



is the density of k states, whether in a cubic or in a 
spherical box |Q [l!|. Equation (||) is fully compatible 
with the development in Refs. Q and JtJ, as it should. 

Since the dipole matrix element tends to a constant in 
the limit /% — > 0, it is easy to see from ^) that the Wigner 
threshold law holds; namely, that lim,j ^o To/wo is finite, 
and nonzero except for an unlikely accident. As a matter 
of fact, our argument about the k — > limit of the dipole 
matrix element was nothing but a recital of a standard 
argument for the Wigner threshold law. Nonetheless, it 
will furnish a relevant piece of the puzzle when we are to 
discuss atom-molecule field theory below. To this end, 
we note from Eqs. (||) and (^|) the equality 
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B. Many atoms in momentum representation 

1. Basic Hamiltonian 

In Ref. H we introduced a phenomenological second- 
quantized Hamiltonian for photoassociation of bosonic 
atoms to (obviously) bosonic molecules, 
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E 



47^ 



S Hk 2 

y + 277 



a k ak 



E 

kk'q 



dkk' • E 



Ah 



J k+ fe'+q 



ak.cik.' + H.c. 



.(12) 



Here m stands for the mass of the atom, m = 2/i. The 
operators a k and 6 k are boson annihilation operators for 
atoms and molecules in the plane wave mode k. The 
Hamiltonian trivially includes the kinetic energy of the 
atoms and molecules. Here photoassociation takes place 
with an absorption (as opposed to induced emission) of 
a photon. By momentum conservation, atoms with wave 
vectors k and k' plus a photon with wave vector q must 
then make a molecule with wave vector k + k' + q. This 
explains the form of the cubic operator product; annihi- 
late the atoms and a photon, create the molecule. 

As is usual in quantum optics, we use a classical 
field to represent the photons. Specifically, the positive- 
frequency part of the electric field reads 



E+(r) 



2 Z-^l 

q 



E r 



(13) 



Both E + (r) and the coefficients E q may in principle be 
slowly varying functions of time, but the leading time de- 
pendence of the electric field oc e~ lu}t has been absorbed 
into the detuning So in a transformation to a rotating 
frame. The seemingly unexpected factor | in the detun- 
ing term correspond to the fact that upon photodissoci- 
ation one molecule produces two atoms, both of which 
generally take away kinetic energy with them. The sign 
of the detuning is chosen in such a way that So > corre- 
sponds to tuning of the laser by the energy TiSo above the 
photodissociation threshold. A quick way to verify this is 
to consider the potential resonance when the (quasi) en- 
ergy (in the rotating frame) for a system of one molecule 
and zero atoms would be the same as the energy for a sys- 
tem with zero molecules and two atoms, much like in the 
Appendix |A|. The total energy of the atoms must equal 
the energy of the molecule plus hSo, which is compatible 
with Eq. (|12|). 

We assume that the optical transition responsible for 
photoassociation is a dipole transition. In the dipole ap- 
proximation, the electronic photodissociation and pho- 
toassociation transitions in a molecule do not (cannot!) 
depend on the propagation direction of light, hence the 
dipole matrix elements d kk / do not depend on pho- 
ton momenta. Besides, by translational invariance, the 
dipole matrix elements must be functions of the differ- 
ence k — k' only. Because of the exchange symmetry of 
the boson operators, the matrix element d kk * may be 
chosen to be symmetric in the exchange of the momen- 
tum indices, hence, an even function of k — k'. And 
finally, since we consider s-wave photoassociation only, 
the matrix element is a function of |k — k'|. 

To pin down the matrix elements d kk /, in Ref. || we 
took the limit of a dilute thermal gas. The known QC 
results are recovered if one expresses the thus far unde- 
fined coefficients d kk / in terms of the matrix elements (^) 
as 



Ikk' 



V2d(i(k-k')) 



(14) 



The y/2 is a consequence of the Bose-Einstein statistics, 
and the factor i follows from the way that the relative 
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momentum must be defined to make it the conjugate of 
the conventional relative position. We have recently no- 
ticed that there is a subtlety associated with this identi- 
fication having to do with the the statistics of the atoms. 
We elaborate in App endix |A| , but meanwhile continue 
according to Eq. (|l4|). 

In sum, if one may treat the atoms and the molecules 
as bosons in their own right, then both the form of, and 
even the numerical coefficients in, the photoassociation 
Hamiltonian are unambiguously determined by simple 
physical considerations. Admittedly, we do not know of 
any conclusive ab initio argument for the bosonic nature 
of atoms and molecules in photoassociation. But neither 
do we know of any for alkali atoms, which nonetheless 
seem to behave like good bosons in current BEC experi- 
ments. 



2. Two-mode approximation 

We now revisit the situation that was the focus of 
Rcf. & an infinite, homogeneous condensate and one 
plane wave of light with photon momentum Tiq. By con- 
servation of momentum, the molecules one may create 
by photoassociation all have the wave vector equal to q. 
On the other hand, as far as momentum conservation is 
concerned, a molecule with the wave vector q may pho- 
todissociate into two atoms, neither of which has zero 
momentum. We will discuss such "rogue photodissocia- 
tion" in more detail below, Sec. VA, and so far ignore 
it. 

All told, we only retain the two modes of the model 
with the annihilation operators a = ao and b 
Hamiltonian reads 



fe q . The 
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(15) 



In the "approximate" equality we have added a constant 
of the motion to the Hamiltonian, a step that has no 
effect on the ensuing dynamics. The QC Rabi frequency 
reads 



lim ' 
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(16) 



where we have used Eq. (||). Without any loss of gen- 
erality, we have chosen k to be real and nonnegative. 
Seemingly alarmingly, the Rabi frequency still depends 
on the quantization volume. We will return to this point 
in Sec. HA. At present, our aim is just to set up the 



second-quantized Hamiltonian in the two-mode approxi- 
mation. The task is completed by noting that 



kg 2 

Am 



(17) 



is the detuning corrected for the photon recoil energy of 
the molecule. From now on we will keep track of this dis- 
tinction, so that S always includes the appropriate recoil. 
In fact, that was already implicitly the case in Eq. (0). 



C. Many-atom field theory 

Our phenomenological Hamiltonian ( |l2| ) was written 
down originally in momentum representation. Neverthe- 
less, as in Refs. || and jl3], it is fairly straightforward to 
convert it into position representation, i.e., into a quan- 
tum field theory. 

Let us introduce atomic and molecular fields in a quan- 
tization volume V. A priori, this volume does not have 
to be the cavity of radius R, as in Sec. II A. Often it is 
actually more convenient to use a cubic box with peri- 
odic boundary conditions. On the other hand, when one 
deals with two atomic fields, it is often expedient to write 
the integrals in the theory in terms of center-of-mass and 
relative coordinates. If the integral over the relative co- 
ordinates cuts off because the integrand tends to zero at 
large distances, it is immaterial what volume is used, as 
long as it is large enough. Thus, when convenient, in 
such relative-coordinate integrals we may still imagine 
the spherical potential well. The short of the story is 
that the quantization volume V refers to the geometry 
that is expedient in the particular context. 

We write the atomic and molecular fields in terms of 
plane wave states as 



e ik ' r a fe , 



v k 



(18) 

The Hamiltonian (|l2j) may be cast as a Hamiltonian den- 
sity of a field theory for these fields. Using the standard 
continuum limit 



V 



(2tt) 3 



d 3 kf(k). 



(19) 



Eq. (|13|), and the properties of Fourier integrals, the part 
of the Hamiltonian density depending on the dipole in- 
teraction becomes 



•V f (r) /d 3 r>(r+ir')d(r')</>(r-ir') + . 



H(r) E+(r ) 

h ~ 2h 

(20) 

The dipole kernel of photoassociation d(r) is given in 
terms of the dipole matrix element of Eq. (H) as 



d(r) 
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It should be noted that our argument contains a subtle 
trick, which is exposed in Appendix |b| 

To make further progress, we study the kernel using 
the qualitative model for the dipole matrix element (Q). 
Some more juggling with Fourier transforms gives 



d(r) 



V2d 



d(a+r)(a+r)i[>*(a + r)-9(a-r)(a-r)if)*(a -r) 
r 

(22) 



where 9 is the Heaviside unit step function. This result 
is only as good as the assumptions of Eq. (0). In par- 
ticular, as Eq. (0) is not valid for large k, Eq. ( p2| ) is 



dubious at short distances. Also, as we already pointed 
out with Eq. (Q), even for small k the model is quantita- 
tively accurate only if the vibrational wave function 4>(r) 
resides outside the range of the atom-atom interactions. 
In spite of these caveats, Eq. (22) gives an exceedingly 
plausible idea of the range Ar of the dipole interaction 
kernel; Ar is of the order of the larger of the spatial ex- 
tent of the vibrational state ?/>(r) and the absolute value 
of the scattering length a. 

The relevant length scale for the atomic field is larger 
than Ar if the energies of all of the atoms relevant for the 
time evolution are smaller than ~ h 2 /m(Ar) 2 . Hence- 
forth we assume this to be the case. Then we may sim- 
plify the field-theory form of the integral as follows, 
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^ f (r)0(r)0(r) + ... . 



(23) 



The second equality follows from the assumedly slow de- 
pendence of the atomic field </>(r) on position r compared 
to d(r), the third from Eq. (Ell) and the properties of 
Fourier transformations, and in the last equality we have 
substituted Eq. (|5|). Comparison with Eq. ( p"T| ) then 
gives the contact interaction form of the photoassocia- 
tion dipole interaction, 



H(t) 



-V(r)ip\r) cj>{r)(j)(r) 



with 



p(r) = ei a rg [d.E + ( r) ] Hm 2nh 2 T Q (r) 



(24) 



(25) 



Here To(r) is the photodissociation rate as per the pre- 
vailing light intensity at position r, and the exponential 
records the phase of the dipole interaction term. 

For reference, we summarize the contact-interaction 
version of the Hamiltonian for photoassociation in terms 
of atomic and molecular fields. While at that, by adding 
a suitable multiple of the conserved particle number com- 
pletely analogously to the operation that we did in the 
chain of equations (15), we move the detuning term to 
the molecules. The result is 
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III. COHERENCE IN TWO-MODE MODEL 

We consider the two-mode model, whose Hamiltonian 
we reproduce here from Eq. (|l5|), 



H 



a) a — ^(b^aa + ba) a)) . 



This is the age-old Hamiltonian for second-harmonic gen- 
eration, with atoms replacing the fundamental-frequency 
light and molecules the second harmonic. Needless to say, 
the literature on the model is extensive. We will cite a 
few relevant examples below. 

Let us take the conserved particle number a^a + 2b% 
to have the value N . The Hamiltonian ( |l5| ) may then 
be restricted to the space spanned by the states \n) = 
\n) m\N —2n) a with n = 0, . . . , N/2 molecules and N—2n 
atoms, and is tridiagonal in that basis. It is easy to find 
both the eigenstates of the Hamiltonian and the time evo- 
lution of any specified initial state numerically. We have 
done so using inverse iteration, and a variation of the 
Crank-Nicholson method. While the problem considered 
in Ref. Q is different, the numerical methods described 
therein work equally well here. 

In this way we have firstly ^l| found the fraction 
of atoms converted into molecules, given that the system 
starts out at time t — with everything in atoms and is 
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driven by a resonant laser, 6 = 0. We found a nonlinear 
analog of Rabi flopping, the system oscillating between 
atoms and molecules. Given that we have a quantum 
system where the spacing between the successive energy 
eigenstates is not constant, it is not a surprise that the 
oscillations collapse, and even revive . 

Interesting as these features are as a matter of princi- 
ple, the drawback remains that Rabi oscillations of oc- 
cupation probabilities are generally not robust, not even 
in quantum optics or laser spectroscopy. We foresee lit- 
tle experimental utility for nonlinear Rabi oscillations. 
Our main message of the analysis of Rabi oscillations in 
Rcf. Q rather is that, for N ^s> 1, the characteristic fre- 
quency scale of the system is not the QC Rabi frequency 
k but 



1.0 



= \/N k= lim 

v ->0 



'2iV7r7rr 
p 2 Vv 



lim - 



* 2itTi Top 

(1 2 Vq 



(27) 

The y/N is evidently a Bose enhancement factor. The im- 
plicit quantization volume has disappeared from the re- 
sult. Instead we have the density of atoms if all molecules 
were dissociated, 



P : 



N 
V 



(28) 



which is a true physical parameter for the system. 

Bose enhancement highlights the role of statistics in 
our analysis. Suppose we were to deal with a nondegen- 
erate gas. Then the occupation probability of each QC 
state is much less than unity by definition, and there is 
no Bose enhancement. The frequency of Rabi flopping 
between a given QC state and a bound molecular states 
would just be cx k cx V~ 1 / 2 , and vanishes in the limit of 
large volume. It is not possible to have Rabi flopping in 
photoassociation of an infinite and homogeneous nondc- 
generate gas, even in principle Jl5| . 

The photoassociation rate for a pair of atoms starting 
in a given QC state is cx n 2 cx 1/V, and also vanishes 
in the limit of a large volume. The saving grace || 0| 
is that, to keep the density constant, the atom number 
N and at the same time the number of candidate atoms 
to photoassociate with any given atom tends to infinity. 
Adding the probabilities for photoassociation due to all 
colliders, the total rate of photoassociation for any given 
atom is cx Nk 2 cx p. This remains a constant in the 
continuum limit, which leads to a finite photoassociation 
rate cx p. 

In contrast, in a BEC the atom number already ends 
up in the transition matrix element, giving the effective 
Rabi frequency Q = y/~N k cx y/p. One does not add rates 
cx k 2 due to different atoms, but all the condensate atoms 
act as a single state that has a finite matrix element for 
photoassociation even in the limit of infinite quantization 
volume. This ultimately facilitates coherence effects like 
Rabi flopping. 

The Hamiltonian (]l^) acts in a rotating frame, and 
there is no manifest physical significance to its eigenval- 
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FIG. 1: Fraction of atoms forming molecules, /, in the 
ground state of the atom-molecule condensate as a function 
of the detuning of the laser from photoassociation threshold, 
S. The long-dashed and short-dashed lines are from the two- 
mode quantum treatment for the respective atom numbers 
N = 10 and N = 100, the solid line is from the semiclassical 
approach embodied in the Gross-Pitaevskii equations. 



ues. Nonetheless, let us call them energies, and the state 
with the lowest energy the ground state. Qualitatively, 
when the laser is tuned far below the photodissociation 
threshold, 6 <C — fi, the a^a term in the Hamiltonian be- 
comes costly in energy. The ground state of the Hamil- 
tonian should then tend to have the system mostly in 
molecules. Conversely, for S 3> fl, the ground state favors 
atoms. We confirm these surmises explicitly in Fig. [II by 
plotting the fraction of atoms converted into molecules, 
the expectation value 



/ = 



N 



(29) 



for the numerically obtained ground state. The difference 
between the dashed curves is the atom number, N = 10 
and N = 100. 

The idea is near that, when the laser is swept from 
a large above-threshold detuning to a large (in absolute 
value) below-threshold value, the system will follow adi- 
abatically, and an atomic condensate is converted to a 
molecular condensates M. Moreover, if Q were the rele- 
vant frequency scale, adiabaticity means that the detun- 
ing must change by in a time of the order or longer 
than ft . In our model we assume that the detuning is 
swept as 



5(t) 



-£,n 2 t. 



(30) 



We expect adiabatic following when £ ~ 1. 

This piece of intuition is correct. In Fig. ^ we repro- 
duce the relevant figure from Ref. Sj. We fix N = 100. 
At an initial time giving S = — 20 f2, we start the sys- 
tem in its ground state, almost purely atoms, and inte- 
grate the time dependent Schrodinger equation numeri- 
cally while sweeping the detuning at two rates £ = 1 and 
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o . fl- 




FIG. 2: Fraction of atoms converted to molecules, / = 2(6^6), 
as a function of instantaneous detuning, S(t), when the de- 
tuning is swept linearly as a function of time according to 
5(t) = — £,Q, 2 t. The figure is for the atom number N = 100, 
for the two sweep rates £ = 0.1 and £ = 1. 



£ = 0.1. The plot gives the fraction of atoms converted 
to molecules as a function of the instantaneous detuning. 
For £ = 1 we expect to be somewhere at the borderline 
of adiabaticity, and actually find a conversion efficiency 
of 0.8. When the laser tuning is swept ten times more 
slowly, £ = 0.1, the conversion efficiency has reached 0.97. 

We believe that we have discovered a feasible method 
for preparing a molecular condensate: start with an 
atomic condensate and a photoassociating laser, then 
sweep the frequency of the laser in a proper way, and 
obtain a molecular BEC. In the general manner of adia- 
batic methods, this "rapid adiabatic passage" should be 
robust. 

Besides its potential utility, rapid adiabatic passage 
is another example of coherent phenomena that occur 
in a condensate but not in a nondegenerate gas. Once 
more, in a nondegenerate gas the frequency scale for 
adiabaticity is k, and the time scale is proportional to 
k" 1 oc V 1 / 2 . Even if nothing else went wrong, in the 
limit of a large sample the time scale for adiabatic atom- 
molecule conversion would tend to infinity. Once more, 
Bose enhancement saves the day by turning the volume 
dependence of the relevant frequency scale into a density 
dependence Ea|. 



IV. FIELD THEORY FOR ALL MODES 

If only one spatial mode has to be considered for both 
the atomic and the molecular condensate, the two-mode 
Hamiltonian ([l5|) is all there is to it. Nonetheless, even 
if somehow the infinite homogeneous condensate could 
be approximated in practice, for whatever reason there 
would always be atoms and molecules present with mo- 
menta that are not included in the two-mode picture. As 
the full Hamiltonian ( [l2| ) mixes momenta nonlinearly, the 
possibility of instability arises and should be investigated. 



We will study deviations from the two-mode picture, 
or the possibility that the (infinite) atomic and molecular 
condensates are not spatially homogeneous. Here we use 
the quantum field version of our photoassociation theory. 
But solving nonlinear quantum field theories is generally 
not a simple matter, and one must approximate. We will 
resort to classical field theory, the same way one proceeds 
when an alkali condensate is analyzed using the Gross- 
Pitaevskii equation. 



A. Gross-Pitaevskii equations 

To derive the Gross-Pitaevskii equations (GPE) for 
this problem, in the usual way we first amend the Hamil- 
tonian by adding to it a multiple of the conserved particle 
number 

N = $3(a£ok+2&£&k) = I d?r [0 t (r)0(r)+2^ t (r)VXr)] . 

k J 

(31) 

Instead of the Hamiltonian density, we then use the 
"Kamiltonian" density 

£(r) = H(r) - ntf (r)</>(r) + 2^(r)^(r)} (32) 

in our calculations. This substitution has no effect on 
the dynamics. At this point the real scalar /i is arbitrary. 
Though this piece of knowledge has no bearing on our 
analysis, in the end the constant [i will be analogous to 
chemical potential in thermodynamics. Hence, \i is called 
the chemical potential. 

Using the standard commutators for boson fields, the 
Heisenberg equation of motion for the molecular field be- 
comes 



iip(r) = [ip(r), 

nv 2 

4m 



^(r) - ±Z>(r)0(r)0(r) . 



(33) 



The transformation to the classical field theory is effected 
by positing that the fields in the equations of motion are 
no longer boson fields, but commuting c-number fields. 
The interpretation is that </>(r) and ij}{r) are the macro- 
scopic wave functions for atomic and molecular conden- 
sates. 

From this point onward we again take the driving field 
to be a simple plane wave with the positive frequency 
part |Ee* q ' r . Second, we scale the atomic and molecular 
fields with the square root of density, yjN/V. Third, 
we incorporate the spatial variation of the electric field 
into the definition of the molecular field. Altogether, we 
define new atomic and molecular fields $ and \& as 



(34) 
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The net results are four. First, the normalization of the 
fields now reads 



- J d 3 r(|d>| 2 + 2|*| 2 ) = l 



(35) 



Second, the coefficient T> gets multiplied by the square 
root of the atom density that would prevail if all 
molecules were to dissociate. Third, all explicit posi- 
tion dependence disappears from the equations of mo- 
tion. The GPE for the rescaled atomic and molecular 
wave functions are 



i<i>(r) 
i*(r) 



hv 2 



2m 

.A- 

4w 

\e iv Sl [$(r)] 2 



$(r) 



h 

2 m 



" i¥ T2$*(r)*(r),(36) 
*(r) 

(37) 



X7 + S-2 



None other than our characteristic frequency scale for 
photoassociation, £1 of Eq. ( p7| ) , is now explicitly the fre- 
quency scale in the field equations as well. The phase 
factor e lcp accounts for the phases of the electric field 
and dipole moment, and will soon prove inconsequential. 
Fourth, a photon recoil term got added to the kinetic 
energy of the molecules. 

While still using full dimensional units, we pause to 
discuss the mathematical symmetries of the GPE. First, 
a trivial phase change of on e of the fields, e.g., <& — > 
e - l ¥>/ 2 <j) ; converts Eqs. (IV A) to the same equations, ex- 
cept that the phase factors vanish from atom-molecule 
interaction terms. Therefore, we drop the phase in the 
interaction term. Second, the GPE have a global gauge 
invariance. If the fields W are a solution, then so are 
the fields §e lv , ^>e 2llp for an arbitrary fixed phase <p. In 
particular, putting ip = 2n, one may see that the equa- 
tions are invariant under the change of the sign of the 
field <&. Third, as a time dependent generalization of the 
gauge invariance, the GPE is invariant under the replace- 
ment 



fi — > fi + Aju . 



(38) 
(39) 
(40) 



Fourth, the GPE is Galilei invariant, in that the replace- 
ments 



*(r,i) -> e lkr $ Ir - — t,t ] , 
m 



#(r,i) -> e^ k r * ( r- — i,t ) , 

m 



^k 2 

2m 
Tik 



(41) 
(42) 

(43) 
(44) 



convert a solution into another solution that corresponds 
to an added momentum Kk. per atom. The nontrivial 



part of the transformation, Eq. ( ff4[ ) is that, due to the 
Doppler shift, the effective detuning changes depending 
on the overall motion of the at om-m olecule system. 

Finally, we express the GPE ( IV A ) in a specific system 
of units, 



to 



1 

fi' 



ro 



h 

mfi ' 



mo = m , 



(45) 



for time, length and mass respectively. Technically, we 
should rename the scaled variables and parameter; say, 
t = tto, 8 = 8 /to- However, we eschew such a heavy 
notation, and continue to use, e.g., the symbol 8 for what 
now stands for a dimensionless number and should be 
properly denoted by 8 = 6/Q. The result is 



[-|V 2 + iiq- X7 + S-2/I] 4- 



2^ 



(46) 
(47) 



The GPE equations are as stated in Ref. @. The 
plan now is to analyze them. Of course, equations of 
the type (IV A) are standard fare in studies of second- 
harmonic generation. The literature once more is exten- 
sive, and we will only give a few specific pointers as we 
go along. 



B. Steady state of GPE 

Atomic and molecular fields <& and \& represent a sta- 
tionary state, one in which the physics does not change 
with time, if and only if their time evolution is solely in 
global (position independent) phase factors, possibly dif- 
ferent ones for <I> and 4". Now consider fields of the form 
$(r, t) = $(r)e~ iwt and *(r, t) = *(r)e" 2iwt for any real 
ui. One sees right away that for this type of evolution, at 
least the exponential time dependence properly matches 
on both sides of Eqs. ( I V A ) . Although we have not been 
able to prove it mathematically, we conjectu re tha t, as- 
suming time independent parameters in Eqs. ( I V A ) , this 
type of time dependence is also the only possible case in 
which the physics is independent of time. 

But then, by virtue of the transformation ( [V A ) , by 



readjusting the chemical potential one can always remove 
the time dependence of the fields entirely in any station- 
ary state. The chemical potential started its life as an 
arbitrary parameter. From now on we make use of the 
arbitrariness and choose the value in such a way that in 
steady state the atomic and molecular fields are literally 
constants in time. 

We shall see shortly that, for any time independent de- 
tuning 8, one may pick a sui table value for the chemical 
potential /i so that the GPE ( IV A ) have solutions <i>o, 
that are constants in both space and time. But by virtue 
of the Galilean transformation ( IVA ), we may then con- 
struct from <E>0j a stationary solution corresponding 
to an arbitrary overall flow of atoms and molecules, so- 
lutions of the form $ oc e zk r and 4 oc e 2lk r . Conversely, 
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we believe that all spatially homogeneous stationary so- 
lutions are such Galilean boosts of the once-and-for-all 
constant solutions $o, V^o- As far as it comes to spatially 
homogeneous steady states, we therefore may, and will, 
without any restriction on generality consider fields that 
are constants in space and time. 

In the context of second-harmonic generation it is well 
known that the GPE may have spatiall y in homogeneous 
solitary- wave type stationary solutions [|23| . By applying 
the Galilean transformation, one may find traveling soli- 
tary waves as well. Here we will make no effort to find, 
let alone classify, solitary- wave solutions to our GPE, but 
proceed as if the homogeneous steady states were all there 
is to it. 

Within the scope of the present paper, the stationary 
solutions thus satisfy 



(S - 2/i)* 



2 



1*0 



2|*n 



(48) 
(49) 
(50) 



where the final equation originates from the normaliza- 
tion (35). Moreover, by virtue of the global gauge in- 
variance, one may choose $o real. Next, because of the 
invariance of the GPE with respect to the sign change of 
$, we may always choose $o to be nonnegative. Further- 
more, if $o > 0, then (|4|) shows that W (like /i) must be 
real. On the other hand, if (f>o = 0, comes with an ar- 
bitrary phase, and may be chosen real. All told, we onl y 
need to find the real solutions <& , *&o, M to Eqs. ( [ V B ) , 
and besides only solutions with $o > need be retained. 

There are three distinct solutions. The trivial one 
reads 



for 6 < a/2 we have 



(51) 



a/g - s 2 - sVeTs 2 



(i 



_ -6-V6 + P 



G 



and for 6 > —a/2 we find 



(52) 



*0 = 



-5 + VWTP 



g 



(53) 



At least two steady states are found for each S, and 
three in the interval —a/2 < 5 < a/2. The ques- 
tion is, which one represents the desired physical steady 
state. Here we attempt to mimic the ground state of the 



quantum- mechanical two-mode model (|15|), and choose 
the stationary solution so that mM 



S < -a/2 : $ = $g, \P = Vq, M = M" 
6>-V2 : $ = $o,* = ^,/i = ^- 



(54) 



The success is evident in Fig. Ill, where we plot side by 



side the fraction of atoms converted to molecules from 
the quantum-mechanical two-mode model and the cor- 
responding GPE approximation / = 2^/ 2 (solid line) as 
a function of detuning. By comparing with the N = 10 
and N = 100 quantum results, it may be seen that the 
agreement between the GPE and the quantum approach 
gets better as the number of atoms is increased. Even 
for an atom number as small as 100 and around the non- 
analytic point S — —a/2 of the GPE approximation (p34|), 
the difference is only on the order of one per cent. 



We conclude by noting that Eqs. (51) and ([52]) together 



give a stationary solution that is the exact mirror image 
of our choice ( |54| ) with the substitution 5 — > —5. This 
corresponds to the state with maximum energy for the 
two-mode quantum system. It is a stationary state just 
as is the minimum, and may be used for atom-molecule 
conversion. The difference is that, in the case of the 
maximum, the detuning would be swept in the opposite 
direction, from below to above the dissociation thresh- 
old. Otherwise rapid adiabatic passage should work es- 
sentially as before. 



C. Small excitations of the system 

So far, while analyzing the GPE, we have achieved 
nothing more than in our studies of the two-mode model; 
rather less, because in our classical field theory we lose 
quantum fluctuations. Nevertheless, we now have the 
tools to analyze the stability of the steady state, and see 
how spatial inhomogeneities evolve in time. 

As in Ref . |l4| , we linearize the GPE around a station- 
ary solution, and then attempt to find eigenmodes for 
small deviations from the stationary case. Both of these 
steps are achieved at once by inserting the Ansatz 



$(r,t) = $q 



e 



(p-r — tut) 



-i(p-r— u)*t) 



(55) 



into the GPE, and only retaining the first-order terms in 
the "small" coefficients w$, v$. Of course, the field \t 
is treated similarly. We need to mix plane waves p and 
— p because the GPE mixes fields and their complex con- 
jugates. With our Ansatz we also retain the possibility 
that the evolution frequency of an eigenmode could be 
complex. 

The Ansatz succeeds if the as of yet unknown evolution 
frequency lu satisfies the eigenvalue equations 

[-5P 2 + A* + w] + ipov<s> + $0^* = 0, (56) 

[-^P 2 + M - w ] v -f + "00 + $o«* = 0, (57) 

[-|(P 2 + 2 P ' q) + 2 /" - $ + + $ u* = 0, (58) 

[~l(P 2 - 2p • q) + 2/z - <5 - w| + $ w* = 0. (59) 
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The characteristic equation is fourth order in w, so in 
principle the solutions can always be written down an- 
alytically. However, here we use Mathematica [pi[ to 
produce results directly numerically, and occasionally to 
extract analytical forms for special and limiting cases. 

The dispersion relations of the excitation modes, lo = 
u)(p), depend on the relative propagation directions of 
the excitation and of light, and also on the size of the 
photon recoil kick. We encompass these dependences into 
a dimensionless parameter, which in terms of the original 
dimensional parameters reads 



V 



' ft(q-p) 2 
mfl p ■ p 



(60) 



Given the parameter r\ and the absolute value p of the 
(dimensionless) propagation ve ctor p , the terms that de- 
pend on photon recoil in Eqs. ( |IV C| ) are replaced as fol- 
lows, 



i(p 2 ±2p-q)^ l(p 2 ±2r,p). 



(61) 



For fixed values of the parameters 5, rj and p, there 
are four excitation modes with four (in general) different 
evolution frequencies lo. If a positive imaginary part is 
encountered in any one of the four frequencies, the corre- 
sponding mode grows exponentially and the steady state 
is unstable. 

Let us take the wave vector characterizing the small- 
excitation mode to be perpendicular to the propagation 
direction of light, r\ — 0. We begin withp = and assume 
| 'So 1 2 < \ i a nontrivial mix of atoms and molecules. We 
find the solutions to the characteristic equation for the 
eigenvalue problem ( IV C| ) 



«?(0) 
lo 2 2 (0) 



2*o 



(62) 
(63) 



When the wave number of the excitation moves away 
from p — 0, of the four frequencies the two evolving con- 
tinuously from ±102(0) remain real. Their dispersion re- 
lations for small p are of the type 0J2(p) — ^2(0) + 2^P 2 , 
with m* > 0, so that these excitations are akin to optical 
phonons. On the other hand, the remaining two excita- 
tion frequencies are either purely real or purely imagi- 
nary, and for small enough p they are always imaginary. 
Such modes do not propagate at all, but grow or shrink 
in place exponentially. 

In fact, in Fig. || which is transcribed from Ref. [fTi) , 
we plot the largest imaginary part among the four fre- 
quencies lo, \^s[loi]\ 7 as a function of the detuning 5 and 
the wave number p of the excitation mode. As above, 
we set 77 = 0. It may be seen that, for any detuning 
5 > — -\/2, an unstable excitation mode is always found. 
The largest imaginary part of an evolution frequency, 
i.e, the largest growth rate of an instability, is encoun- 
tered at 5 = -0.154496 and p = ±0.771324, and equals 
3?(w) = 0.24256. 



. 25 R 
0.2 ' 

0.1 
0.05 




5 



-2 



FIG. 3: Largest imaginary part among the four evolution 
frequencies of small perturbations of the joint atom-molecule 
condensate is plotted as a function of the laser detuning 5 
and the wave number of the plane-wave like perturbation p. 
In this figure the excitation mode propagates perpendicularly 
to the propagation direction of light, so that £ = 0. 



The analogous instability is naturally known in second- 
harmonic generation ; and goes under the rubric "modu- 
lational instability" [£5], g6| . 

We next come to the effect of the direction of the wave 
vector of the mode on instability. In Fig. [| we plot the 
largest positive imaginary part of a mode frequency found 
for any p as a function of the parameter n. The curves 
are labeled with their corresponding fixed detunings S. 
In our studies of this kind, the largest growth rate of the 
instability was always found for r\ = 0. For a given de- 
tuning, the small excitations whose momentum direction 
is perpendicular to light propagation always present the 
most unstable scenario. 

Finally consider the largest growth rate of the instabil- 
ity (among all p) as a function of the detuning. As may 
be inferred from Fig. |^, with increasing detuning <5 ~ 1 it 
decreases and occurs at smaller momenta, i.e., at larger 
size scales. In fact, for i5> 1, the maximum growth rate 



12 



0.25 




0.00 



FIG. 4: Largest imaginary part in a small-excitation evolu- 
tion frequency encountered for any value of the wave number 
of the excitation p is plotted against the parameter £ that 
characterizes the component of the photon recoil momentum 
in the direction of the propagation vector of the small exci- 
tation p. The explicit definition of £ is given in Eq. (po[). 
Different curves are for different fixed detunings S, as identi- 
fied in the legend. 



of an instability R(S) — max{SJ[u;(p, 5)]} p behaves as 



R(8) 



26 463 +0( -S^ 



(64) 



and the corresponding position of the maximum pu (8) , 
the value of the momentum such that Q[u>{p, 5)] < 
3[w(pm, 5)] for all p, goes like 



Pm{5) 



1 



47 
16 <5 4 



8 6 



(65) 



Of course, as the largest growth rate of the instability 
seems to occur for r\ = 0, this was our choice in Eqs. (|64|) 
and (H). 

In sum, we have found that for 8 < the steady 

state of the atom-molecule system is stable, and for any 
8 > — \/2 it is unstable. Although we have reported only 
on a specific stationary solution (|54|), possibly one out of 
three, the modes we have not considered explicitly are 
all unstable. But 8 = — V2 is also the watershed, in 
that below 8 = — v2 the steady state is all molecules 
(<f>o = 0), and above 8 = —\J~2 some atoms are involved. 
All told, the steady state with everything in molecules is 
stable for 8 < —y/2, but any steady state involving atoms 
is always unstable. The most unstable situation occurs 
around 8 ~ 0, with about an equal mixture of atoms and 
molecules. For 5 ^> 1 atoms take over, and the time scale 
for the instability grows longer. 



D. Fate of unstable system 

Linearized stability analysis has shown that a joint 
atom-molecule condensate is unstable in the presence 
of photoassociating light, i.e., there are small-excitation 
modes that grow exponentially. But when a small devi- 
ation from steady state grows exponentially, eventually 
it is not a small deviation anymore and the linearized 
analysis breaks down. To investigate the fate of the sys- 
tem once the instability has set in, we integrate the full 
GPE numerically as in Ref. Jl4| , Our aim is a qualita- 
tive demonstration, so we proceed in 1 + 1 dimensions, 
one spatial coordinate x and time t. However, there is 
nothing in our method that would not immediately work 
in higher spatial dimensions. The restriction is merely a 
matter of computer time. 

Specifically, in our algorithm we discretize a stretch L 
of the line into equidistant points Xi, and seek to repre- 
sent the fields at these discrete points only. For conve- 
nience, we use periodic boundary conditions, so that the 
value of all functions repeats over L. The GPE is inte- 
grated over a time step h in two moves. First we ignore 
the position derivatives in the GPE altogether. This im- 
plies that the stripped-down version of the GPE is local; 
for each Xi, $(xi) and ^>(xi) only couple to $(xi) and 
^>(xi). We integrate the local version of the GPE over 
the time h separately for each Xi as two coupled differ- 
ential equations using a second-order Runge-Kutta step. 
To prevent a numerical drift of the norm, we complete 
the initial step of the algorithm by normalizing the fields 
analogously to Eq. ((3^). Second, we ignore anything but 
the position derivatives in the GPE. The resulting par- 
tial version of the GPE is nonlocal, but in exchange it is 
linear and does not mix the fields <!> and ^ . To integrate 
over the same (sic!) time step as in the initial part of 
the algorithm, we first take the Fourier transform of the 
result from the first part using the Fast Fourier Trans- 
formation (FFT). In Fourier space position derivatives 
become local, so, to propagate the fields over the step h, 
we simply multiply their Fourier transforms by what is 
now the local, exact (within FFT), linear time evolution 
operator. Finally, in preparation for the next step, we 
transform back to real space. 

This split-step algorithm is an obvious variation of the 
time honored split-operator methods for linear partial 
differential equations p7j , and has been described be- 
fore at least by the group of Firth p8[ . Nonetheless, it 
comes with a fair dose of heuristics. It is therefore grat- 
ifying that we have been able to verify a good rate of 
convergence by using successively smaller time steps to 
integrate over a fixed interval of time. 

We present an example of our simulations in Fig. 
showing the absolute square of the atomic field |$| 2 as 
a function of position x and time t. We use 128 points 
Xi. Because of the periodic boundary conditions, the left 
and right edges of x wrap around and are actually the 
same. The plot is for the parameters 8 = 0, T) = 0, the 
range of position x is 24.6282 (in units of ro) correspond- 



13 




around, and tend to join when they collide. As far as 
we can tell, within the present model only one big, dense 
clump will remain in the end. 

One might wonder what is the mechanism behind the 
clumping. We present a heuristic guess. We are not talk- 
ing about a thermodynamic system, so minimization of 
energy is a dubious principle to begin with; and besides, 
we are dealing with quasienergies in a rotating frame, not 
real energies. With these warnings out of the way, let us 
surmise that the system nevertheless attempts to mini- 
mize its energy. The atomic and molecular condensates 
make something akin to the two-level system in quantum 
optics. When light is added, one gets a dressed two-level 
fi system. The energy of the lower one of the two dressed 
4 states decreases with increasing Rabi frequency, which is 
1 the product of the electric field strength of the laser and 
o the dipole moment. But the analog of Rabi frequency 
for the two-level system of atomic and molecular conden- 
sates, i7, is also proportional to the square root of density, 
so the present system may also decrease its energy by in- 
creasing the density. Maybe this is what the instability 
is about. 



EXPERIMENTAL CONSIDERATIONS 




FIG. 5: An example to demonstrate the modulational insta- 
bility of the joint atom-molecule condensate. Atom density 
|<I>| 2 is plotted as a function of the spatial coordinate x and 
time t, given that the system starts out at time t = in the 
steady state plus a small amount of random noise. 



ing to three wavelengths of the most unstable excitation 
mode for these parameters, and time t runs from to 127 
(in units oi to). The unit of |$| 2 , atom density, is such 
that for a homogeneous gas with everything in atoms, 
the density would be |$| 2 = 1. The system starts at time 
t = in the steady state appropriate for 8 = 0, rj = 0, 
except that we add a small amount of Gaussian noise to 
each of the points specifying the initial state. Figure || 
is otherwise the same as Fig. 2 in Ref. except that 
a different seed was used for the random number gener- 
ator that added the noise. As befits an instability, this 
innocuous change has lead to a totally different quanti- 
tative result at long times. 

From the numerical simulations we see that the nature 
of the instability is such that the atoms and the molecules 
together combine into dense clumps. These clumps move 



The models we have discussed until now have been 
rather rudimentary. We next take up two types of com- 
plications that may come up in experiments. First, in 



Sec. VA enters an angle that is in principle included 



in our many-body Hamiltonian, although we have not 
yet considered it expressly: photodissociation of conden- 
sate molecules to states outside of the atomic condensate. 
It turns out to limit the rate at which one can achieve 
coherent atom- mole cule conversion in adiabatic passage. 
Second, in Sec. V B we discuss a number of aspects that 
have so far been missing from our models: spontaneous 
emission from the photoassociated state, atom-atom in- 
teractions, trapping of atoms and molecules, and various 
level shifts. Spontaneous emission can be ameliorated by 
resorting to a two-color photoassociation scheme. Pro- 
vided that photoassociation is speedy enough, which we 
believe is technically possible, the rest of these complica- 
tions may be minor nuisances rather than dominant fea- 
tures of an experiment. What it takes to ma ke ph otoas- 
sociation speedy enough is the subject of Sec. |V Q where 
we discuss the characteristic Rabi frequency for phot oas- 
sociation for various alkalis. Finally, in Sec. 
analyze the published experiment of Ref. JlF 
viewpoint of coherent photoassociation. 



VD 



frc 



we 



A. Rogue photodissociation 

As we noted already, in the case of photoassociation 
of an infinite homogeneous condensate, momentum con- 
servation uniquely determines the state of the ensuing 
molecule. The converse, however, does not hold. When 
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a molecule photodissociates into two atoms, momentum 
conservation does not force the atoms to return to the 
atomic condensate. Bose enhancement favors recombi- 
nation of atoms with the condensate; the characteristic 
frequency is Cl for both photoassociation and photodis- 
sociation between atomic and molecular condensates. 
Nonetheless, atoms winding up elsewhere are lost for co- 
herent photoassociation. We have coined the term "rogue 
photodissociation" for photodissociation processes that 
send atoms outside the atomic condensate. 

One might think that energy conservation in the cy- 
cle of photoassociation and photodissociation is the ad- 
ditional constraint that guarantees that the atoms return 
to the condensate. But this need not be a compelling ar- 
gument. Any time dependence in the system interferes 
with energy conservation. For instance, suppose that 
photodissociation proceeds to the noncondensate states 
at the same rate Tq (per atom) as it would in the case 
of a nondegenerate gas of molecules, so that after a time 
~ Tq 1 coherent photoassociation ceases. The photodis- 
sociation rate Tq in itself furnishes a time scale such that 
energy has to be conserved only to within KTq. The time 
evolution involved in nonlinear Rabi flopping would also 
interfere with energy conservation. 

We present here a rudimentary model for rogue pho- 
todissociation for the special case when the detuning is 
swept in order to convert an atomic condensate to a 
molecular condensate. The key assumption is that we 
may employ the standard Markov approximation in the 
analysis of photodissociation. This entails that rogue 
photodissociation has no memory, but is characterized 
at each instant of time by a rate of exponential decay. 
Such an assumption seems dubious in particular when 
the laser is tuned to the close vicinity of the photodis- 
sociation threshold [ p9| . However, we know of no near- 
threshold case of this kind in which the breakdown of the 
Markov approximation has proven relevant in an experi- 
ment. 

Evidently, only the condensate mode, one of very many 
atomic modes, is strongly affected by Bose enhancement. 
We take rogue photodissociation to proceed at the rate 
that would be appropriate for a nondegenerate gas at 
the given detuning. Second, we model the dependence of 
the photodissociation rate on detuning using the Wigner 
threshold law, so that we write 



(66) 



For convenience we have chosen the photoassociation fre- 
quency scale Cl as the reference detuning for photodis- 
sociation rate; Tq is the photodissociation rate for the 
detuning S — CI. As the third quantitative element of 
the model, we take the probability that a given atom is 
in the molecular condensate to be twice the square of 
the molecular field amplitude as solved from the classical 
field theory, and normalized as in (|5(]). Explicitly, in di- 
mensional units and for 8 > 0, this probability is found 



from Eq. (153^ as 



Pm(S) = - 




(67) 



Suppose now that the detuning is swept as 8 = —£Cl 2 t, 
as in our rapid adiabatic passage example. Ignoring the 
depletion of the condensates due to the very same rogue 
photodissociation, we find the total probability for rogue 
photodissociation 



P = 



<5(t)>0 



dtP M [8{t)]vJ 5 -^- 



Cl 



£Cl 



2 3/4 

3V4 



dr 



(68) 



where the numerical constant has the value a — 4.03197. 
Except for the numerical factor, this is the same expres- 
sion we already used in a qualitative estimate in Ref. . 

Since the photodissociation rate grows linearly with 
light intensity, cx I, and the photoassociation character- 
istic frequency CI is proportional to the field strength of 
the laser, cx y/l, in the end rogue photodissociation wins 
out as the light intensity is increased. Qualitatively, when 
P = 1, rogue photodissociation has overtaken coherent 
conversion of atoms to molecules. To study this border- 
line case, we first set P = 1 in Eq. J68| ) and solve Tnas 
a function of CI. We then insert the result into Eq. ( |27|) , 
thus eliminating To. Moreover, the velocity vo in Eq. fl27[ ) 
is then the relative velocity of the dissociated atoms cor- 
responding to the detuning of the laser that gave the 
photoassociation rate Tq, in this case 8 = Ct. Therefore 
we have fg/2/i = hCl, and vq may be eliminated as well. 
We finally solve for the borderline value CI as a function 
of the problem parameters. After simple manipulations 
the ensuing characteristic frequency scale for photoasso- 
ciation may be written 



CI = 



2/3 



(pX 3 ) 2/3 e R . 



(69) 



Here we have introduced A, wavelength of the light di- 
vided by 27r, and the familiar recoil frequency for laser 
cooling 



h 



2mX l 



(70) 



The main finding is that rogue photodissociation re- 
stricts the light intensity that one may profitably use for 
adiabatic atom-molecule conversion. This means that 
there is also a maximum usable photoassociation fre- 
quency, or a minimum possible time scale for adiabatic 
atom-molecule conversion. The way we have written our 
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estimate ( |69| ) , the frequency scale is provided by the pho- 
ton recoil frequency, and the corresponding time scale is 
in the ballpark of e^ 1 . The density dependence of pho- 
toassociation is encapsulated in the parameter pX 3 , the 
usual dimensionless parameter that governs the coupling 
of light with matter in a dense medium. For present-day 
condensates pA 3 ~ 1 is a reasonable rule of thumb. Fi- 
nally, we have a numerical constant that depends on the 
rate of sweeping of the detuning, but which may also be 
set equal to one in a rough estimate. Altogether, when 
in need of a qualitative number for the photoassociation 
frequency fi, we resort to Q ~ 6r. 

Although our estimate of the minimum time scale for 
coherent photoassociation was developed for rapid adia- 
batic passage, we believe that (with £ ~ 1) it also applies 
to Rabi flopping. This is because the dimensional pa- 
rameters of the problem are the same in both cases. In 
fact, as it comes to dimensional quantities, the minimum 
time scale for coherent photoassociation is equivalcntly 
written fcj 



This is the essentially unique quantity with the dimension 
of time that can be put together using the quantities 
characterizing a homogeneous, noninteracting, quantum 
mechanical, zero-temperature BEC; density, atom mass, 
and fi. 



B. Physics missing from model 

1. Spontaneous emission 

We have discussed a one-color model for photoasso- 
ciation. The physical drawback is that, where there 
is a strong dipole matrix element for photoassocia- 
tion/dissociation driven by external light, there is also 
a strong dipole matrix element for spontaneous emission. 
For instance, if a photon is absorbed in photoassocia- 
tion, then there is also a reverse spontaneous decay of the 
molecule. The molecule may end up in bound vibrational 
states, either in the same electronic manifold where the 
atoms started from, or in some other electronic manifold. 
Alternatively, the photoassociated molecule may decay 
back to a dissociation continuum in a process known as 
radiative escape. Either way, usually the probability is 
small that the system returns to the same two-atom state 
in which is started. After each process of spontaneous 
emission, two atoms are typically lost for any profitable 
use. 

The analogous problem of an unstable excited state 
is standard fare in quantum optics, and so is the solu- 
tion: add another laser-driven transition from the spon- 
taneously decaying state to a stable state. In the same 
way, two-color Raman photoassociation, a free-bound 
transition followed by a bound-bound transition of the 
molecule, may take place between (nearly) non-decaying 
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FIG. 6: Schematic three-level scheme for photoassociation. 
The levels a, b and g stand for atoms, primary photoassociated 
molecules, and final stable molecules, respectively. The figure 
is drawn in the rotating frame, so that the Rabi frequencies 
k and x representing the laser fields are independent of time, 
and the intermediate and two-photon detunings A and 8 are 
drawn as appropriate level energies. 

atomic and molecular states fj], ^2|, [53|. Here we study 
Raman photoassociation as a means of achieving an effec- 
tive two-mode scheme; genuine three-mode phenomena 
such as STIRAP are discussed elsewhere |D| 

Our present notation for this scheme is sketched in 
Fig. ^. Suppose the first step of photoassociation takes 
place with absorption of a photon, then one sets up an- 
other laser to force, say, induced emission from the pri- 
mary photoassociated state to a (more) stable bound 
molecular state. We call the Rabi frequency in the second 
step Xj the detuning of the laser from resonance in the 
second transition A, and the spontaneous decay rate of 
the primary photoassociation state T s . In this context S 
stands for the two-photon detuning, the total energy mis- 
match for light-induced transition from the initial atoms 
to the final stable molecular state, including appropriate 
photon recoil corrections. 

Let us model our three-level A scheme using a variation 
of the two-mode model (|l5|) as 

^ = (5.9+3 + (A+8)tfb- x{b ] g+g^b) - \n{tfaa+ba)a)) , 

(72) 

where g, for "ground" , is the annihilation operator for 
the stable molecular state. Once more, we have added a 
multiple of the conserved particle number, this time in 
such a way that the atoms have zero energy. The term 
oc x describes transitions between the bound molecular 
states. Spontaneous losses from the intermediate state 
are ignored for the time being. 

The Heisenberg equations of motion from Hamilto- 
nian (|72|) are 

g = -i5g + ixb, (73) 
b = -i(A + S)b + ixg + \ima, (74) 
a = ina/b. (75) 

We next eliminate the intermediate state adiabatically 
with the assumption that A is the largest evolution fre- 
quency in the system. We thus formally set b — 0, and 
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obtain 



2. Interactions between atoms and molecules 



Inserting this into Eqs. (|73j) and (|75|), we have 



i . «x 

t:1 — aa, 
2 A 
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(76) 



(77) 



(78) 



These are Heisenberg equations of motion for the effective 
Hamiltonian 
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g*aa) 



(79) 

The effective Hamiltonian describes a two-level system 
with the two-photon detuning 5 and two-photon Rabi 
frequency kx/ A in lieu of the usual detuning and Rabi 
frequency. There are two additional twists to the story. 
First, the two-photon resonance experiences a light shift 
— x 2 /A, an old acquaintance from quantum optics. Sec- 
ond, we have an effective atom-atom interaction propor- 
tional to k 2 /A. An analogous interaction in a nondc- 
generate gas was discussed earlier in Ref. pjj| . Other 
than these tweaks, everything we have said about the 
two-mode model applies as before. 

The adiabatic elimination (76) has a dark side hidden 
by a notational trick, namely, that it does not preserve 
boson commutators. Had we written the right-hand side 
in Eq. ( |75| ) as ba^ instead of a 1 b and substituted (|76|), the 
atom-atom interaction term in the effective Hamiltonian 
would have displayed the operator ordering aaaJa) in- 
stead of aJafaa. The difference is immaterial for a large 
atom number, the limit we are studying anyway, but a 
more careful investigation of the adiabatic elimination 
would be in order if the atom number were not large. 

As excessive care with operator products is not war- 
ranted, we write from the adiabatic assumption the num- 
ber of atoms in the intermediate state b qualitatively as 



g 1 g- 



1 2A") 



rjLV( a t a ) 2 



(80) 



Let us scale the operators by y/N, i.e., write b = \/N0, 
and so forth. Then, within a factor of two, the quantum 
expectation value {j3'0) is the probability that an atom 
is in the intermediate state, and so on. Equation J8C| ) 
becomes 



(4 



V7+(^) («M 2 . (81) 

The fraction of atoms lost per unit time to spon- 
taneous emission from the intermediate state equals 
2r s (/?t/3).The intermediate detuning A suppresses losses 
by a factor oc 1/A 2 , whereas the effective Rabi frequency 
scales as 1/A. In principle, and at the present level of 
the physical model, it is possible to get rid of the harmful 
spontaneous emission to any desired degree by increasing 
the intermediate detuning. 



Atoms interact among themselves, molecules interact 
with molecules, and atoms even interact with molecules. 
For a dilute gas, the atom-atom interaction is often de- 
scribed by the effective two-body potential 



f/(ri,r 2 ) 



4nh 2 c 



S(ri - r 2 ) , 



(82) 



where a is the s-wave scattering length for the atoms 
as before. One may write analogous models for atom- 
molecule and molecule-molecule interactions. If atom- 
atom interactions were suspected to be a factor, one could 
add to the field theory the usual two-body atom-atom 
interactions as 



Ti-AA — 



2nh 2 a 



0t(r)0t( r )0( r )0( r ) 



(83) 



and so on. Inelastic collisions, such as quenching of the 
molecules by collisions, may also prove important. At a 
phenomenological level, they could be described by using 
a complex scattering length for atom-molecule collisions. 

Nonetheless, we have considered neither elastic nor in- 
elastic collisions explicitly in this paper. The motivation 
is mainly pragmatic. Photoassociation is the novelty of 
this work anyway. Second, at this point in time virtually 
nothing is known about the scattering lengths for cases 
other than atom-atom interactions. Third, as pointed 



out in Sec. VA, we anticipate a characteristic frequency 
scale for photoassociation, f2, to be of the order of pho- 
ton recoil frequency 6r of laser cooling, say, ten kilohertz. 
This is larger than a typical frequency scale associated 
with collisions, Aithap / m, in many of the present alkali 
experiments. Photoassociation should dominate the ac- 
tion at least over short time scales. 

More formally, we see from Eqs. (|69|) and (71) that the 



maximum usable photoassociation frequency scales with 
atom density as p 2 ^ 3 , while the rate of binary collisions 
scales as p. In principle and at this level of modeling, it 
is always possible to make photoassociation win out by 
decreasing the density. 

Of course, not all of our discussions are for short times 
only. Notably, in the case of the instability of a joint 
atom-molecule condensate, it may well happen that col- 
lisional interactions eventually play a role in the clump- 
ing. We plan to return to collisional effects in a future 
publication, inasmuch as something worthwhile emerges 
from this front. 



3. Trapping of atoms 

In the current alkali vapor experiments one does not 
see infinite homogeneous condensates, but the conden- 
sate is ordinarily confined to a magnetic trap with a 
(practically quadratic) potential V(r). This may be 
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taken into account in the field theory by adding a term 
in the Hamiltonian density, 

HvA(r) = ^(r)V(r)<f>(r). (84) 

A similar additional term would describe trapping of 
molecules. 

In fact, the beauty of the field theoretical formulation 
is that it does not depend on any given one-particle ba- 
sis to describe the motion of the atoms and molecules. 
The photoassociation term was originally discussed us- 
ing plane-wave states, which makes the derivation easy, 
but at the level of field theory there is no manifest vestige 
of plane waves anymore. Even if we add trapping poten- 
tials, there is no need to tamper with the photoassocia- 
tion term. This should be contrasted with the more del- 
icate situation that emerges if one tries to consider pho- 
toassociation directly using the eigenstates of the trap, 
or indeed some states that would take into account both 
trapping and atom-atom interactions. 

Once more, if the photoassociation frequency scale is 
of the order it is still vastly larger than the typi- 
cal frequency scales associated with magnetic trapping 
of atoms, and the corresponding length scale for pho- 
toassociation, A, is far smaller than the size of a typical 
condensate. Over short times the condensates behave lo- 
cally as if they were homogeneous. One just applies the 
theory of an infinite condensates at each local density, 
and averages the results over the trap. 

On the other hand, there are cases in our formulation 
where the trapping will matter. For Q ~ en the modu- 
lational instability may well set in as in a homogeneous 
condensate, but the motion of the atoms and molecules 
due to the trapping forces will certainly have a long-term 
effect on the atom-molecule clumps. We do not discuss 
this issue here, but plan to return to trapping in a future 
publication. 



4- Level shifts 

We have already mentioned a few mechanisms that can 
change the position of the photoassociation resonance. 
Atom-atom and molecule-molecule interactions alter the 
energy per atom or per molecule, and thereby modify 
the resonance condition for photoassociation. Since the 
atom-molecule ratio conversely depends on the detuning, 
the makings of bistability and hysteresis are in princi- 
ple there. We have also discussed the light shift and the 
many-body shift in a two-color, three-mode configura- 
tion. 

One more shift we have brought up before |(| but 
not yet in this paper, arises because the dissociation con- 
tinuum is not flat. Given the initial state, the dipole ma- 
trix element (or more precisely, the square of the dipole 
matrix element per unit energy) depends on the final 
continuum state. The result is that the photodissocia- 
tion rate picks up an imaginary part. That, of course, 



amounts to a shift of the photodissociating state with re- 
spect to the continuum. The shift is proportional to light 
intensity, and in a qualitative estimate is comparable to 
the photodissociation rate; see, e.g., Refs. j34|. 

Moreover, there are additional shifts due to the pres- 
ence of each and every discrete state dipole-coupled to 
the initial state. The sum of all light shifts is actually fi- 
nite only because the dipole coupling eventually tends to 
zero when one goes high enough in the energy of the cou- 
pled states. The implication is that, a priori, all dipole 
coupled states, even those off resonance by perhaps sev- 
eral photon energies, have to be considered explicitly. If 
one finds a significant contribution to the light shift from 
one far-off resonance state, chances are that one has to 
consider all of them. 

We will not attempt to address continuum and non- 
resonant light shifts explicitly. Nonetheless, on the ba- 
sis of the atomic case discussed in Refs. p3|, we believe 
that, at least above the photodissociation threshold, the 
light shift should be reasonably independent of where ex- 
actly the laser is tuned. This would mean that, in the 
rapid adiabatic passage type atom-molecule conversion, 
the light shift merely gives a constant bias to the detun- 
ing, and is virtually inconsequential. 

C. Numerical examples 

1. From rate to Rabi frequency 

One-color photoassociation has been analyzed by vari- 
ous groups |35|] , in particular for alkalis at finite tempera- 
ture. A typical outcome is the photoassociation rate 1Z v i 
(in s _1 ) or the photoabsorption rate coefficient a^, A (in 
cm 5 ) for a given bound level v' of the excited electronic 
state. The latter quantity is independent of two experi- 
mental parameters, namely atom density p and laser in- 
tensity /. The rate of photoassociation TZ v i is obtained 
via 

= p<pa% A , (85) 

where the photon flux (photons/scm 2 ) is given by ip — 
I /fkOL, wl being the photon frequency. 

But we also know from our earlier work j^, |?J that 
the photoassociation rate in a thermal sample is given in 
terms of our detuning 8, temperature T, and photoasso- 
ciation rate F as 

1l v , = p\ 3 D F , (86) 

where 

- Get) (87) 

is the usual thermal de Broglie wavelength, albeit cal- 
culated using the reduced mass of the colliding atoms. 
In what follows, we write the detuning parameter for 
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photoassociation as S = Woo — ujt, — A„< , where hu^ is 
the asymptotic energy difference between the electronic 
curves and A„/ is the red-detuning of level v' from its 
asymptote. The binding energy of the molecular state is 
thus equal to hA v >. Combining Eqs. J85| ) and (^) with 
Eq. (27), we have an expression for the characteristic fre- 
quency of coherent photoassociation to the level v', 



V V V 2X3 D 

In all of our discussion below we choose the detuning 
in such a way that HS = ^ksT, which gives the cor- 
responding resonance velocity v = yjh^Tj\L. Strictly 
speaking, Eq. (27) requires the limit v — > 0, or equiva- 
lently, 5 — > 0. However, we always assume, without ex- 
plicitly checking this assumption, that the temperature 
is already low enough to bring the system into the region 
of validity of the Wigner threshold law. The quotient 
T/v in Eq. ( ]27| ) has then supposedly reached the v — > 
limit. 

We have already introduced the usual density param- 
eter for light-matter coupling A 3 p to characterize atom 
density, and the recoil frequency en as the frequency 
scale. In the same vein, we write the Rabi frequency 
for photoassociation as 



(89) 



where I v > is the characteristic light intensity for coher- 
ent photoassociation. In explicit numbers, we find from 
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(90) 



which also displays the units used to express atomic mass 
m, temperature T, wavelength of light A, and photoab- 
sorption coefficient a^, A . 



2. Calculated photoassociation rates 

One can estimate the photoabsorption rate coefficient 
for a pair of atoms. At low temperatures, only s-wave 
scattering contributes to the process, and one finds 0, 



^KT)^— A^e-^|^(M)| 2 , (91) 

In the low-temperature limit, the dipole matrix clement 
|-D„/(M)| 2 can be approximated by 

\D V ,(E)\ 2 = \{v',J=l\D\hS,J=Q) , 

" @)lA)| 2 (a-i^')X< (92) 



TABLE I: Photoabsorption rate coefficients a p , A (in cm 5 ) 
and characteristic intensities /„/ (in W cm -2 ) for levels v' 
with the best simultaneous Franck- Condon factors with the 
continuum and the highest-lying level of the lower electronic 
state, and for levels corresponding to detunings near 1 cm -1 . 
The vibrational number v' as well as the corresponding bind- 
ing energy A„/ (in cm - ) are also given. Calculations were 
performed for T=l mK. 
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with E — h 2 k 2 /2/i = h5, where Dq = |d| is the asymp- 
totic dipole moment, a is the scattering length, R v i is the 
classical outer turning point of the excited level v', and 
S v ' is a dimensionless parameter representing the frac- 
tion of the bound wave function contained in the last 
node §§§. 

For example, given 7 Li atoms in a triplet state at 1 
mK, a detailed calculation j37j showed that the excited 
level v 1 — 58 has the best Franck-Condon factor with 
the continuum and the highest-lying bound level v" — 
10 of the lower triplet electronic state, with af/i 58 = 

2.0 x 10~ 32 cm 5 pSfl . The corresponding characteristic 
intensity is I V '=5g = 33Wcm -2 . 

Beyond this, in Fig. ^| we give the calculated rate co- 
efficients for high-lying vibrational states for both stable 
isotopes of Li, for both the singlet and the triplet excited 
states. Correspondingly, Table | presents a few numer- 
ical examples of the rate coefficients and characteristic 
intensities. 

Much larger photoassociation rates and correspond- 
ingly smaller characteristic intensities can be obtained 
if the excited levels v' are closer to the dissociation 
limit (i.e., have smaller binding energies). In fact, the 
v' = 58 level is deeply bound, with a binding energy of 
hA v , =58 = 5.5 x 10" 4 hartree (or 122.22 cm" 1 ||). For 
example, v' — 86 with a binding energy of 0.87 cm 
has a rate coefficient 1050 times larger, namely af/i 86 = 

2.1 x 10 -29 cm 5 , which gives the characteristic intensity 
I v '=86 = 0.03 2 Wcm~ 2 . The high-lying excited levels v' 
have larger rate coefficients because the overlap of the 
excited bound wave function and the continuum ground 
wave function is larger. For high v', the overlap scales 
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|). Notice that 

1/2 



80 60 

a v < (cm 1 ) 

FIG. 7: Variation of photoassociation rate coefficient a^A for 
high- lying vibrational states. In (a), the r Li singlet transi- 
tions over a range of 600 cm" 1 exhibit oscillations mapping 
the nodal structure of the continuum wave function, with val- 
ues varying over five orders of magnitude. Plots (b) and (c) 
show a shorter range of detunings for singlet and triplet tran- 
sitions, respectively. The singlet scattering lengths of both 
isotopes are positive and a^A decays rapidly from its large 
values near A ~ cm -1 , by nearly three orders of magnitude. 
The triplet scattering lengths are both negative, and a^A de- 
creases slowly for small detunings, until the first nodes of the 
continuum wave functions are reached. 



This equation is similar to Eqs.flSq 

since Ad oc the rate scales as 1/T and A^, 1 ^ 2 for 

small detuning. For lithium at 0.140 mK (correspond- 
ing to the Doppler temperature), assuming a laser inten- 
sity of 1000 W/cm 2 (so that k = 0.8 x 10 10 s -1 : see 
pj"| for details), and a detuning A v > ~ 1 cm" 1 ~ 30 
GHz (or 4.5 x 10 -6 hartree), we have i?„< = 135 ao and 
X D = 1489 a = 7.875 x 10" 6 cm. For a density p = 10 11 
cm" 3 , this gives p\ z D = 4.88 x 10" 5 and, neglecting the 
scattering length a, one gets 1Z V > ~ 9.6 x 10 4 s" 1 . If we 
scale the rate to T = 1 mK, we get 1.3 x 10 4 s -1 , less 
than twice the value 7140 s" 1 obtained from Eq. ( |85| ) 
using the calculated value of ar^/£ 86 given above. No- 
tice that many assumptions are made in these estimates: 
nonetheless, expression ( |94"| ) gives good order of magni- 
tude for the photoassociation rate for small detuning and 
low temperatures. 

Using the Doppler temperature as a typical temper- 
ature, and averaging the rate using a linewidth corre- 
sponding to 5 MHz [[ll] , one finds typical rates for small 
detunings by scaling the numbers of Table || with the 
appropriate temperature, density, laser intensity, and de- 
tunings, according to 



lo po 1 V A v > 



(95) 



Here, J = 1000 W/cm 2 , p = 10 11 cm" 3 , A = 1 cm" 1 , 
and To and TZ are listed in Table |l[ Correspondingly, 
a PA scales as 



like A" 1/2 . 

Similar results for sodium and cesium are available. 
For small binding energies, the rate coefficient a p A at 
1 mK is of the order of 10~ 28 cm 5 for Na, and 10" 27 
cm 5 for Cs |Q. These translate to characteristic inten- 
sities as small as /iWcm" 2 . Notice that for lithium, the 
rate coefficients for singlet transitions are smaller than 
for triplet transitions, reflecting the sign of the scatter- 
ing length: this is especially significant for 6 Li, where the 
triplet scattering length is negative and enormous |lq| . 

Other expressions based on semi-classical treatment 
have been developed jit], |il| . For example, in jy], the 
rate for s-wave contribution is 



o — 

p\ A D e "bt 



3SA V . 



sm 2 [k(R v , - A{k))}n 2 ; (93) 



PA -^o / A PA 

a *' r"V a7 



(96) 



The characteristic intensities should then scale with the 
square root of the detuning A v / , so that we have 



Iv> = 



A,,/ 



Ar 



■To- 



(97) 



The values of A, (p and a PA are listed in Table || 
To compare the rates for the various alkali metals, it 



TABLE II: Approximate scaled photoabsorption rates TZ and 
rate coefficients a PA for weakly bound molecular states. The 
values of TZ, and T are from Pillet et al. Q, I = 1000 W 
cm -2 , po = 10 11 cm ~ 3 , and Ao = 1 cm" 1 . 



where n = d ■ E/2h is the Rabi frequency of the laser of 
intensity /, A(k) represents a shift from the free solution, 
and A v i (again) is the detuning of the level v' from the 
dissociation limit. This expression is valid for s-wave 
scattering and at low detunings. In the limit k — > 0, 
A(kf{) — > a, and the maximum of 1Z V / is found to be at 
hS ~ k B T/2, so that 



TZ v i ~ pX 



D<- 



a) 2 . (94) 



Atom T TZ A ip 

mK 10 4 s" 1 nm 10 21 cm" 2 nn 



PA 
5 



Li 0.140 45 

Na 0.240 22 

K 0.140 25 

Rb 0.140 13 

Cs 0.125 10 



671 
589 
766 
780 
852 



3.4 
3.0 
3.9 
3.9 
4.3 



1.3 x 10" 27 
7.3 x 10" 28 
6.6 x 10" 28 
3.3 x 10" 28 
2.3 x 10" 28 
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TABLE III: Approximate photoassociation rates 1Z, rate co- 
efficients a PA , and characteristic intensities I v , for the alkalis. 
The fixed parameters are Jo = fOOO W/cm 2 , pa — fO 11 cm -3 , 
Ao = f cm -1 , and T=f00 pK,. Also listed are the reference 
densities A -3 (in cm~ 3 ) and frequencies en (in 27rkHz). 



Atom 


n v , 


a p , A 

V 






A" 3 






10 4 s" 


1 cm 5 




mWcm -2 


cm" 3 


27rkHz 


7 Li 


63 


1.8 X 10" 


27 


3.7 


8.21 x 10 14 


63.3 


23 Na 


53 


1.7 x 10" 


27 


0.47 


1.21 x 10 15 


25.0 


39 K 


35 


9.2 x 10" 


28 


0.18 


5.52 x 10 14 


8.72 


87 Rb 


IS 


4.6 x 10" 


28 


0.069 


5.23 x 10 14 


3.77 


133 Cs 


12 


2.9 x 10" 


28 


0.039 


4.01 x 10 14 


2.07 



is convenient to express them for the same parameters. 
Assuming I = 1000 W/cm 2 , p = 10 11 cm" 3 , T = 100 
/iK, and Ao = 1 cm" 1 , we obtain the values listed in 
Table H. The photoassociation rates vary between 1.2 x 
10 5 s-Mpr a PA ~ 2.9 x 10~ 28 cm 5 ) for Cs and 6.3 x 
10 5 s" 1 (or a PA ~ 1.8 x 10~ 27 cm 5 ) for Li. From the 
expressions for 1Z v i , one expects the rates to scale like 
l/2/i. If we were to multiply the rates by M, the mass 
number, we notice that, except for Li, the scaled rates 
indeed are similar. The variations left are due to slightly 
different Rabi frequencies |tj| . 

Finally, one has to be careful when using the values 
of Table [II, in which the effect of scattering lengths are 



not taken into account. These effects can be significant, 
like in the case of 6 Li, 85 Rb or 133 Cs, where large nega- 
tive scattering lengths induce larger rates. Also, at larger 
detunings, one probes deeper region of the excited elec- 
tronic state, and shorter distances of the lower state con- 
tinuum wave function, where the exact nodal structure 
will play an important role (see jl7|, |36[ [57| and Fig. |t]) . 

As already not ed, the values for the saturation inten- 
sity in Table III are to be construed as estimates only. 
Also, as they are, they apply to a very weakly bound 
molecular state (A„' = 1 cm -1 ). Such high- lying states 
are likely to be easily perturbed by atom-molecule colli- 
sions, so that in practice one might resort to more bound 
molecular states. Nonetheless, these saturation intensi- 
ties are remarkably low, occasionally much lower than 
our initial generic estimate of lOWcm -2 ^ Q. From 
the experimental viewpoint this is encouraging, as the 
requirements on laser intensity are greatly moderated. 
The flip side is that the usable intensity appears to be 
limited by rogue photoassociation; for densities such that 
A 3 p ~ 1, the maximum photoassociation Rabi frequency 
is of the order of the recoil frequency . For heavier 
alkalis the saturation intensities may indeed be low, but 
small recoil frequency and slow coherent photoassocia- 
tion are the prices to pay. 



D. Discussion of an experiment 

In the experiment of Wynar et at, Ref. jl6|, the au- 
thors studied photoassociation of a 87 Rb condensate in a 
two-color Raman configuration. From our standpoint, 
the reported data fall into two categories. First, the 
width and shift of the two-photon resonance line was 
studied as a function of the intensities of the two lasers. 
We are going to use these measurements to estimate the 
bound-bound Rabi frequency. Second, the authors mea- 
sured the width and the shift of the same resonance line 
as a function of the density of the condensate. At low in- 
tensities, elastic atom-atom and atom-molecule collisions 
turn out to have a major effect on the resonance parame- 
ters. However, these effects were modeled quantitatively 
in Ref. fl6j] , and one may then get at the quantities per- 
taining to photoassociation. We use this type of data to 
determine the characteristic frequency of coherent pho- 
toassociation. 



1. Adapting the theory 

We first cast our theoretical approach in a form that 
allows for direct juxtapositions with the analysis of ex- 
perimental data in Ref. Jig . W e ig nore elastic collisions 
because a comparison with Ref. 1 16 turns out to be possi- 
ble even without considering them explicitly, and inelas- 
tic collisions because their presence was not conclusively 
established in Ref. |16| ]. 

We begin by adding the spontaneous decay of the in- 
termediate state at the rate V s as a nonhermitian term in 
the Hamiltonian for the A system, Eq. w%. This gives 



? = Sg Jf g+(A + Spb- 
n 



2 



(98) 



It is known in quantum optics that introducing an imagi- 
nary part to the energy of the decaying state gives correct 
results as long as the unstable state decays irreversibly 
only to states that are outside of the state space included 
in the model. 



We then repeat much of the analysis of Sec. VB 1 for 
the Hamiltonian (|98|), with a number of tricks and ap- 
proximations. First, we take the limit of large intermedi- 
ate detuning, |A| ^> T 8 , \S\. Second, we assume that the 
bound-bound Rabi frequency \ is much larger than the 
free-bound Rabi frequency fi. This implies that in the es- 
timates of the line widths and decay rates, such as those 
following from Eq. (|8l|), only the x 2 terms need be kept. 
Third, we ignore light-induced atom-atom interactions. 
Fourth, we scale the the operators as before, a = yNa, 
etc., which brings out the photoassociation frequency O 
explicitly. Fifth, after deriving the equations of motion 
for the atomic operator a and final-state molecule oper- 
ator 7 by adiabatic elimination of the intermediate-state 
operate (3, in the spirit of semiclassical approximation we 
replace the operators a, /3, and 7 with c numbers. At this 
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stage we have, in analogy with Eqs. (|77|) and fl7q), the 
equations 



.(, x 2 \ r sX 2 
1 aJ 7 ~ 2A 2 " 

A ' 



, i • n X > 



<T^> (") 
(100) 



To develop the rate approximation, we first note that 
the probability that an atom belongs to the stable state 
g, P g = 2|7| 2 , satisfies the equation of motion 

P 9 = - T -jfp 9 + ^(C-C*). (101) 
The "coherence" C = a 2 j* has the equation of motion 



C = i^[P g P a -\Pl\- 



2A 2 



C, (102) 



where, in turn P a = \a\ 2 is the probability that an atom 
remains in the system unassociated. 

We could go on and derive a corresponding equation 
of motion for P a . However, here we assume that P a ~ 
1, and that correspondingly P 9 « 1. We solve for the 
coherence C adiabatically by setting (7 = 0, and insert 
the result into Eq. (101) to obtain a rate equation for the 
population of the stable molecular state g, 



Pa 



A 2 9 





2 


rr.x 2 1 


A 




2A 2 





2 

+ 








2A 2 



P z 



(103) 



This displays a shift of the ground state by 



2. Comparison with experiments 

Wynar et al , Ref. jl6| , discuss the light shift of the res- 
onance line. In their case the frequencies of the two lasers 
are quite close so that both of them cause line shifts and 
broadenings, but this is included in their analysis. Trans- 
lating back to our case with only one laser frequency on 
the bound-bound transition, the comparison between our 
formulation and Ref. [fl6| reads 



A 



A! ' 



(106) 



where I2 is the intensity (in W cm' 2 ) of the bound-bound 
light, Ai = A = 2tt x 150 MHz is the intermediate detun- 
ing, and (3 = 3.5 x 10 9 m 2 W -1 s -2 is a parameter that we 



deduce from their Fig. 
formula 



3a. This gives the approximate 



X 



0.94 



Wcm -2 



x 2ttMHz. 



(107) 



Wynar et al, Ref. [fl6| , also present an analysis of 
line shapes at low intensity. The rate coefficient for 
two-photon photoassociation on resonance, Kq = 9 x 
10" 14 cm 3 s _1 , is one of the fitted parameters. The reso- 



nance photoassociation rate from our Eq. ( |103| ) and the 
corresponding expression in Ref. jlqj , ignoring inelastic 
processes and with n = p, boil down to 



O 2 

— = nKn . 



(108) 



Given the spontaneous decay rate of this molecular level, 
r s = 12 x 27rMHz p6| , and noting that their cited value 
of Kq was for the intensity of the photoassociation laser 
I\ = 0.51 Wcm -2 , we find the scaling formula 



a lincwidth 



2 



2A 2 



(104) 



(105) 



and a corresponding direct decay of bound molecules at 
the rate jl- 

With th e as sumption P 9 « 1 we have made the rate 
equation (103) essentially one way, i.e., ignored transi- 
tions from molecules back to the atoms. The semiclassi- 
cal approximation in general requires that there are many 
molecules in a quantum state, so that the molecules make 
a condensate. However, we believe that, by ignoring tran- 
sitions from molec ules back to atoms, we have divorced 
the rate equation (103) from any coherence requirement 
for the molecules. As the experiments do not claim a 
molecular condensate, this is essential for the analysis 
that follows. 



n ~ 5.8 



p 



h 



10 14 cm~ 3 Wcm -2 
the characteristic 



x 2tt kHz . 



(109) 



intensity I v i = 



which gives 
0.08 Wcm" 2 . 

In the experiment, the binding energy of the primary 
photoassociated state was 23 cm -1 . In contrast, for 
At,/ = 1 cm -1 , Table III gives the characteristic intensity 
0.07 mWcm~ 2 . The theoretical and experimental char- 
acteristic intensities differ by three orders of magnitude, 
but they are also for different binding energies. The dif- 
ference in characteristics intensities is no cause for con- 
cern, as we do not know at all how the detuning 23 cm -1 
for Rb is related to the minima of the photoassociation 
rates like those shown in Fig. ^ for Li. 



Reaching coherence 



Suppose that the density of the 87 Rb gas equals p = 
~ 3 , a factor of two larger than quoted in Ref. JIB], the 
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laser intensities are h = I2 = 10Wcm~ 2 , typical val- 
ues in photoassociation experiments, and the interme- 
diate detuning is as in Ref. |^6|. Then the two-photon 
Rabi frequency becomes fix/ A = 1 x 27rkHz, while the 
effective linewidth is x 2 T s /2A 2 = 2 x 27rkHz. For co- 
herent phenomena we need a Rabi frequency at least 
comparable to the linewidth, a condition that is not 
quite satisfied. But if one were simply to take the in- 
tensities 1\ — 100 W cm" 2 for the photoassociating light 
and I2 = 1 Wcm~ 2 for the bound-bound light, Rabi fre- 
quency remains unchanged and linewidth goes down by a 
factor of 100. One is then in principle deep in the regime 
of coherent photoassociation. Also, the effective Rabi 
frequency would be smaller than the recoil frequency, so 
that rogue photoassociation need not yet be a fatal prob- 
lem. 

On the other hand, eyeballing from Fig. 2 of Ref [^6) , 
the linewidth due to elastic atom-atom collisions would 
be of the order of 5 x 27rkHz, so elastic collisions would 
interfere with coherent photoassociation. Rabi oscilla- 
tions would be seriously impeded, rapid adiabatic pas- 
sage maybe less. One could compensate by lowering the 
density by two orders of magnitude and at the same time 
increasing the intensity of the photoassociating laser by 
another two orders of magnitude to keep the effective 
Rabi frequency constant, but with lowering of the den- 
sity one would wind up in trouble with rogue photoasso- 
ciation. 

Experimentally, one can think of varying the intensity 
of the two lasers, the intermediate state, the intermediate 
detuning, and the density of the sample. The conditions 
one must watch out for are that effective Rabi frequency 
be larger than the effective damping rate and line broad- 
ening due to elastic (and inelastic) collisions, yet not so 
large as to cause excessive rogue photoassociation. If the 
parameters could be varied freely, within our physical 
model a solution is always possible. Unfortunately, one 
must live with some practical limitations on, say, laser 
intensities. The present experiment of Ref. jl6| is tanta- 
lizingly close to coherent photoassociation, but to really 
get there may require a careful optimization of the exper- 
imental parameters and a solid understanding of what is 
technically feasible in a given experiment. We will not 
attempt to enter such discussions. 



VI. CONCLUDING REMARKS 

We have developed theory for coherent photoassocia- 
tion of a Bose-Einstein condensate of atoms, which typ- 
ically leads to a Bose-Einstein condensate of molecules. 
Phenomena analogous to coherent optical transients in 
few-level systems, such as Rabi flopping and rapid adia- 
batic passage between atomic and molecular condensates 
are expected. Coherent phenomena depend on the Bose 
enhancement of the dipole moment matrix element, and 
cannot occur in a nondegenerate gas, at least not in the 
thermodynamic limit. Bose-Einstein condensation and 



photoassociation are both currently front-line research 
and the marriage thereof might still be a technological 
challenge, but coherent optical transients in photoassoci- 
ation should eventually be feasible in experiments. 
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APPENDIX A: DIPOLE MATRIX ELEMENT 

In our bare-bones example we consider the situation in 
which we either have one molecule at rest and no atoms, 
or no molecule and two atoms in states k and — k. Then 
the relative momentum becomes fik. This peculiarity, 
and also the reason for the factor i in Eq. (jl4|) , is because 
the relative momentum for two entities with momenta pi 
and p 2 should be defined as ^(pi — P2) to make it the 
conjugate of the usual relative position — r 2 . We thus 
have the possible state vectors 

h£>= (W + E^i^-k) 10} , (Al) 

where [3 and ctk are complex amplitudes to be deter- 
mined, and |0) stands for the vacuum of atoms and 
molecules. The states for and a_k are the same, so 
the sum over k runs only over half of the possible k val- 
ues; say, those with k x > 0. 

Following Ref. B, we write the time dependent 
Schrodinger equation from the Hamiltonian (|l^) in terms 
of the coefficients (3 and ctk, assuming a plane wave of 
light. In order to avoid an inessential complication, we 
ignore photon recoil and simply set q = 0. The result is 

— -5 Q )a k + l ^-P, (A2) 

^ = i E^^. (A3) 

k 

The QC approach discussed in Refs. and JjJ dis- 
plays precisely the same structure, evidently identifying 
d • E/2ft as the QC Rabi frequency k of those papers. 

There is a catch, however. Because of the Bose- 
Einstein statistics, states with the relative momenta for 
two atoms k and — k are the same, so the density of states 
for the relative motion is half of what it was for Maxwell- 
Boltzmann atoms. To obtain the same photodissociation 
rate, we therefore have to make the matrix element a fac- 
tor of y/2 larger than in the case of distinguishable atoms. 
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This is the \/2 in Eq. (|l4|). The rate of photodissocia- 
tion is proportional to the square of the matrix element 
and picks up a factor of two, which compensates for the 
missing half of phase space density. 

The second-quantized Hamiltonian is chosen so that 
it gives the right photodissociation rate. What we did 
not fully realize in Ref. || is that the same factor of \/2 
also leads to twice the photoassociation rate that one 
would obtain for Maxwell-Boltzmann atoms, given the 
Franck-Condon factor from the standard calculations or 
as deduced from the photodissociation rate. It appears 
that in the literature the rate for s-wave photoassocia- 
tion is usually calculated low by a factor of two. Be- 
sides, by a simple extension of the statistics arguments, 
p-wave photoassociation should vanish for bosons. Arc 
these problems? 

We think not. In photoassociation experiments no par- 
ticular care is usually taken to polarize the atoms. They 
are not all in the same internal state, and therefore not 
all photoassociation events are between indistinguishable 
atoms. This should reduce the factor-of-two boson en- 
hancement of the photoassociation rate for even partial 
waves, and call forth odd partial waves. The remaining 
discrepancies with the conventional calculations may well 
be within experimental and theoretical uncertainties. 

It would be of some interest to investigate photoassoci- 
ation of a polarized low-temperature gas experimentally. 
We believe that the statistics effects are real, and should 
be observable. After all, few colleagues seem to have a 
problem with the assertion that there are no s-wave col- 
lisions for polarized fermions. Photoassociation experi- 
ments in a polarized gas should make an interesting test 
of our phenomcnological Hamiltonian. 



APPENDIX B: KERNEL FOR 
PHOTOASSOCIATION IN FIELD THEORY 

The process we followed in the main text in deriving 
the integral kernel d(r) [as in, e.g., Eq. (23)] was a phe- 
nomenological mix between formally pure field theory 
and collisional physics. 

In fact, to obtain the dipole matrix elements dkk' in 
Eq. (n3) , we carry out an integral of the form 



d kk , = d(k - k') = d / d 3 r i,* (r)<^ k _ k , (r) , (Bl) 



where ip is the bound-state molecular wave function and 
</>k-k' describes the state of two atoms with the relative 
momentum h(k— k'). Now, in the transformation to field 
theory (|l8|) we used plane waves as the wave functions 



for st ates with a given wave vector. In the matrix ele- 
ment (Bl) we should correspondingly use plane waves to 
describe the relative motion, as in 



c-k'(r) 



1 



i(k-k')-r 



V 



(B2) 



However, in our derivation we did not not use plane 
waves, but distorted plane waves that take into ac- 
count atom-atom interactions, and only asymptotically 
(at large distances) corresponds to the relative momen- 
tum k — k'. 

This makes a difference. For instance, if the matrix 



element (Bl) is computed using pure plane waves, the 
final photoassociation kernel reads 



d(r) = V2dip*(r), 



(B3) 



instead of something akin to (|22|). Depending on the 
value of the scattering length, there could be a sub- 
stantial difference in the overall numerical factor also in 
the contact-interaction form for atom-molecule coupling, 
Eq. (p^). The question is, which method is "correct"; 
pure or distorted plane waves? 

First of all, atoms do interact, which in general will 
have an effect on photoassociation. One who wishes 
to analyze photoassociation quantitatively will have to 
take atom- atom interactions into account. A field the- 
orist pursuing pure plane waves would have to consider 
atom-atom interactions explicitly, and work out the con- 
sequences in photoassociation ab initio. This is probably 
a major chore. Our hope is that, with our distorted-plane 
wave trickery, we have phenomenologically captured most 
of the effect of atom-atom interactions on photoassocia- 
tion. 

On the other hand, if we use the distorted plane waves 
to calculate the photoassociation matrix element and 
in addition introduce a contact interaction model for 
atoms as in ( |33"| ) , the suspicion arises that we are double- 
counting atom-atom interactions. Further approxima- 
tions, such as the contact-interaction form for photoas- 
sociation or classical field theory, could either exacerbate 
or ameliorate the double-counting. 

Disentangling atom-atom interactions and photoasso- 
ciation might make an interesting theoretical exercise, 
but we do not attempt it in this paper. Instead, we pro- 
ceed with a model that allows for both distorted-plane 
wave matrix elements for photoassociation, and an ex- 
plicit atom-atom interaction. In this way, we at least get 
the limiting cases of photoassociation in a dilute gas and 
atom-atom interactions in the absence of photoassocia- 
tion basically right with little effor 
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